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Abstract 

Plasma instabilities can play a fundamental role in quark-gluon plasma equilibration in the high 
energy (weak coupling) limit. Early simulations of the evolution of plasma instabilities in non- 
abelian gauge theory, performed in one spatial dimension, found behavior qualitatively similar to 
traditional QED plasmas. Later simulations of the fully three-dimensional theory found different 
behavior, unlike traditional QED plasmas. To shed light on the origin of this difference, we study 
the intermediate case of two spatial dimensions. Depending on how the "two-dimensional" theory 
is formulated, we can obtain either behavior. 



I. INTRODUCTION AND RESULTS 



It is a difficult theoretical problem to fully understand how quark-gluon plasmas come to 
local equilibrium in the context of heavy ion collisions. In fact, equilibration is not yet fully 
understood even in the formal limit of weakly-coupled quark-gluon plasmas (equivalently the 
limit of arbitrarily high energy collisions). By investigating the rates and interplay of individ- 
ual scattering processes, Baier, Mueller, Schiff and Son []J attempted to analyze equilibration 
in the weak coupling limit, producing what is known as the bottom- up scenario of quark- 
gluon plasma equilibration. However, interesting collective behavior in the form of plasma 
instabilities turns out to severely complicate the problem 0]. The relevant instabilities are 
known in traditional plasma physics as Weibel or filamentary instabilities @, 0, 0]- Over the 
last few years, there has been a variety of work on the dynamics of plasma instabilities in the 

weakly-coupled limit of non-abelian gauge theory @, B 0, B M, El, (H, H H H H E3- 

Plasma instabilities initially grow exponentially quickly, and they create large color mag- 
netic fields which could speed up local isotropization and thermalization of the initially 
non-equilibrium plasma by scattering plasma particles into random directions. To puzzle 
together the theory of equilibration in the weak coupling limit, it is important to under- 
stand just how large these instabilities and their associated magnetic fields grow. It has 
been of particular interest to understand the similarities and differences between the evo- 
lution of instabilities in non-abelian gauge theory plasmas vs. traditional abelian plasmas 

As we'll briefly review, simulations have been previously devised to test whether non- 
abelian effects would limit the growth of plasma instabilities compared to abelian plasmas. 
An example of the results 0, [ij is shown in Fig. [TJ This figure shows the growth in time 
of the magnetic energy of long-wavelength fields associated with the instability. In these 
simulations, the instability has been seeded with a very small initial amplitude. The dotted 
line shows the abelian case, corresponding to simple exponential growth of the instability 
(with a rate that can be computed in perturbation theory). This energy growth is fueled by 
stealing energy from the particles in the plasma, and it must eventually stop once the fraction 
of energy stolen becomes significant. In these simulations, that scale is beyond the top of the 
graph. (In fact, the back-reaction on the plasma particles is intentionally ignored in these 
simulations, providing an unlimited source of energy, so that one may cleanly disentangle 
whether or not non-abelian effects can limit the instability growth earlier.) 

For computational simplicity, the first simulations of the non-abelian theory were re- 
stricted to one spatial dimension 0] , in a specific sense reviewed later. From rough reasoning 
based on the one-dimensional theory, Ref. 0] had earlier conjectured that non-abelian inter- 
actions would not stop exponential growth of the instability. This conjecture was apparently 
supported by the one-dimensional simulations, an example of which is shown by the dashed 
line in Fig. [TJ The growth starts out looking like the abelian case (which is expected from 
perturbative analysis of the instability). Eventually, there is a deviation when the fields get 
large enough for non-linear self-interactions in the non-abelian theory to become important, 
but the subsequently growth is again exponential. The conjecture proved wrong for the full 
three-dimensional theory, however, shown by the solid line in Fig. [TJ After non-linear inter- 
actions become important, there is a spurt of rapid growth similar to the one-dimensional 
case, but then instability growth stops in three dimensions. There is continued linear growth 
in the magnetic energy at late times, shown in the linear plot of Fig. [Tb, but this is due to 
the unstable modes pumping energy into slightly higher momentum, stable modes rather 



2 




FIG. 1: The magnetic field energy associated with instabilities as a function of time for abelian 
theories (dotted line), non-abelian theories restricted to one spatial dimensions (dashed line), and 
non-abelian theories in three spatial dimensions (solid line). The vertical axis is shown as both 
(a) logarithmic, and (b) linear. These particular simulations are taken from Ref. [g], but similar 
one-dimensional simulations were performed earlier in Ref. 0]. The non- Abelian simulations are 
based on SU(2) gauge theory. 



than by continued growth of the unstable modes [11 . 

The purpose of this paper is to investigate the root cause of the different behavior in one 
and three dimensions. Our tool will be the intermediate case of two dimensions. Before 
explaining our results, however, we need to review the three-dimensional theory being sim- 
ulated and how one defines its reduction to lower dimensions. Following the review, we will 
present our two-dimensional results. We can get either one-dimensional or three-dimensional 
behavior, depending on how the two-dimensional theory is constructed. Then, in Sec. [TTl we 
offer a general interpretation of our results in terms of whether all polarizations of the gauge 
field are gauged in the dimensionally-reduced theory. For the sake of concreteness, we also 
discuss a specific example of this interpretation in the context of Nielsen-Olesen instabilities, 
which could possibly play a role in the fate of non-abelian plasma instabilities. Finally, we 
offer our conclusions in Sec. IHIl Technical details of our simulations are relegated to an 
appendix. 



A. Review 

Weibel instabilities occur when the momentum distribution of a non-equilibrium plasma 
is anisotropic, the expansion and collision rates are small compared to the instability growth 
rate, and there is a clear separation of scale between the momentum of typical plasma par- 
ticles and the wave numbers associated with the resulting instability. The problem can be 
studied by simulating a non-abelian analog of the classical Vlasov equations of traditional 
plasma physics. In the weak coupling limit, one can justify the preceding assumptions and 
the application of these classical equations to the plasma instabilities in the bottom-up pic- 
ture |2j . In the bottom-up picture, the energy stored in the plasma particles is parametrically 
larger than the energy associated with the onset of significant self- interact ions of magnetic 
fields associated with the instability. The Vlasov equations consist of (i) a collisionless 
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Boltzmann equation describing the evolution of the phase space density f(x,p,t) of hard 
quarks and gluons in the presence of a soft gauge field A(x,t), and (ii) Maxwell's equation 
for the soft gauge field with sources given by the hard particles. In order to study how large 
the instabilities grow, it is sufficient to linearize these equations in fluctuations Sf(x,p,t) 
of the hard particles about an initial (homogeneous) distribution fo(p), while retaining the 
non-linear structure of the non-abelian gauge field A. This approximation treats the plasma 
particles as an infinite reservoir of energy for instability growth. Schematically (suppressing 
sums over particles types), the equations are of the form 1 

(A + v ■ D x ) Sf + g (E + v x B) ■ V p / = , (1. la) 

D v = f = gvt R J-0-v»5f. (1.1b) 

Above, Sf is in the adjoint representation, D is the corresponding covariant derivative, tn 
is a group factor, v counts particle spins and types, 2 and = (l,p). For ultra-relativistic 
systems, one can put these equations in a more convenient form by integrating over the mag- 
nitude \p\ of momentum to replace Sf(x,p), which depends on six phase-space coordinates, 
by a distribution W(x, v), which depends on only five since v = p lives on the unit sphere. 
W a (x, v) represents the net adjoint color of all particles moving in direction v at point x. 
The corresponding Vlasov equations are 

(D t + v- D x ) W + m^E ■ (V„ — 2v) — B ■ (v x V v )]Q{v) = 0, (1.2a) 
DvF i*> =f= [ V »W, (1.2b) 



where Q(v) is determined by the angular dependence of the initial particle distribution fo{p), 
normalized so that the angular average of Q(v) is one. The mass scale above is set by 
the amplitude of / as 

m ~ = ^/ (2^^- (L3) 

For moderately anisotropic plasmas, this is the order of magnitude of instability wave- 
numbers, instability growth rates, plasmon masses, and Debye screening, all of which we 
refer to as soft physics. Soft magnetic fields become non-perturbative when B ~ m^/g, 
with corresponding energy density |5 2 ~ m^/g 2 . 

In this paper, we focus on the case of moderately anisotropic, oblate velocity distributions 
which are axi-symmetric about the z axis. In particular, we will study precisely the same 
angular distribution Q(v) studied in Ref. [8j. A full study of the bottom-up scenario will 
require simulating extremely anisotropic distributions 0, 14, Tf|, which is computationally 



more difficult. However, moderately anisotropic distributions will be adequate for the goal 
of this paper to understand the origin of the differences between one and three dimensional 
simulations. For the sake of numerical simplicity, our simulations are all based on SU(2) 
gauge theory rather than SU(3) QCD. We do not have any reason to expect qualitative 
differences between the two. 



1 We work in ( — | — | — (-) metric convention. 

2 For example, for a plasma of gluons, v — 2 and £r = 3. 
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The "one dimensional" version of equations (11.11) and (11.21) consists of assuming that 
fields and distributions vary in space only in one direction z, as A^z^t) and 5f(z,p,t) [or 
W(z, v, t)]. Particle momenta p, however, are still treated three-dimensionally, and all three 
spatial polarizations of the gauge field (A x , A y , A z ) are included. This type of dimensional 
reduction of the full three dimensional theory is known in traditional plasma physics as 
a 1D+3V theory, where the ID signifies that x has been reduced to one dimension, but 
the 3V indicates that velocity (momentum) is still treated three-dimensionally. Fully three- 
dimensional simulations are 3D+3V. Note that if the transverse polarizations A x (z,t) and 
A y (z,t) were not included in one-dimensional simulations, one could not describe magnetic 
physics at all and so could not investigate Weibel instabilities. 

B. Two-Dimensional Results 

Fig. [2] shows our results for two-dimensional (2D+3V) simulations. The solid line repre- 
sents the same kind of dimensional reduction that was done in the one-dimensional case: all 
polarizations are retained, but fields and distributions depend on only two spatial dimen- 
sions: 3 

A fl = A (t (y,z,t), n = x,y,z,t; (1.4a) 
6f = 6f(y,z,p,t), (1.4b) 

or W = W(y, z, v, t). [The labeling of directions with respect to the two-dimensional plane 
is shown in Fig. [3] for reference.] In this case, we find behavior similar to one-dimensional 
simulations: instability growth continues exponentially even after non-abelian interactions 
become important. However, the dashed line shows what happens if we eliminate the po- 
larization A x which lies outside of the two-dimensional plane, by setting A x = in our 
two-dimensional simulations. We now obtain behavior similar to three-dimensional simula- 
tions instead: the energy growth at late times is linear rather than exponential. 




m t m t 



3 Though we have listed /i = t in (|1.4ap for the sake of generality, the simulations are carried out in Ao = 
gauge. We will only discuss gauge-invariant results. 
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FIG. 2: As Fig. [T] but for two-dimensional (2D+3V) simulations. The solid line is for the straight- 
forward dimensional reduction of (jl.4p . The dashed line is a simulation that does not include the 
out-of-plane polarization A x . 




FIG. 3: Axis labeling for two-dimensional (2D+3V) simulations. The spatial plane x is the yz 
plane, shown on the left, but particle momenta p live in three dimensions, and the initial distribu- 
tion Q(v) of particle velocities is axi-symmetric about the z axis, as depicted on the right. 

We can even produce both behaviors in one simulation, shown by the solid line in Fig. 
IU by initializing the seed field for the out-of-plane component A x much smaller than the 
other components. The solid line shows the total magnetic field energy, which, for times 
moot ^ 80, behaves similar to the A x = simulation shown by the dashed line of Fig. [2j 
initial exponential growth followed by linear growth. The dotted line in Fig. H] shows the 
energy associated with A x , which we take to be \B 2 + \B\ = ^(D y A x ) 2 + ^(D Z A X ) 2 in 
two dimensions. The importance of A x continues to grow throughout. When it eventually 
catches up to the other components at m^t ~ 100, the linear behavior changes to the 
continued exponential growth we saw earlier in simulations where A x was initialized on the 
same footing as the other components. 

We have shown simulations that start from tiny initial seed fields for the instability. 
Similar results, shown in the Appendix, apply to two-dimensional simulations that start from 
large, non-perturbative initial configurations: the late-time growth of magnetic energy is 
exponential (like one-dimensional simulations) if the out-of-plane polarization A x is included 
in the simulation, and linear (like three-dimensional simulations) if it is not. 

II. INTERPRETATION 
A. General 

We can now offer a general characterization of what seems to determine whether non- 
abelian effects stop Weibel instability growth in 1D+3V, 2D+3V, and 3D+3V simulations 
(at least for the moderately anisotropic particle distribution we have simulated). Non-abelian 
effects will prevent late-time exponential growth if and only if no gauge-field polarizations 
are included outside of the subset of spatial dimensions simulated. In ID, we have no 
choice; unless such polarizations are included, we cannot simulate magnetic fields and so 
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FIG. 4: As Fig. [2] but for a two-dimensional simulation with very small initial conditions for the 
out-of-plane polarization A x . The solid line is the total magnetic energy |-B 2 ; the dotted line is 
the energy density ^B y + \B\ associated with A x . 

cannot study the Weibel instability at all. In three dimensions, we automatically include all 
polarizations. In two dimensions, we have a choice, with the results described above. 

Formally, the difference between gauge field polarizations in and out of the subset of spa- 
tial dimensions is that those within the subset (like A y and A z in two dimensions) are gauge 
fields in the dimensionally reduced theory, while those outside (like A x ) are instead adjoint- 
charge scalars. For the two-dimensional theory, for instance, the Yang-Mills Lagrangian 
is 

\F«pF™P + \{D a A x ) a {D a A x ) a , (2.1) 

where a and j3 are summed over the 2+1 dimensions (y, z, t). There is a significant difference 
between gauge fields and scalars. Gauge fields can always be locally transformed away so 
that A = at any particular point, and scalars cannot. In our two-dimensional theory, the 
squared amplitude A a x A a x of the two-dimensional scalar A x is locally gauge invariant and can 
have physical consequences. 

As an example, consider the two-dimensional theory, and imagine that a component of 
A has grown very large due to the instability. Let's focus on distance scales small compared 
to the wavelength of the unstable modes. As a first approximation, treat A as constant over 
these small distances and compare and contrast two situations: (i) a large, constant gauge 
component A y or A z , and (ii) a large, constant scalar component A x . In the first case we 
can gauge away to A y = A z = 0, and so there is no effect at this order of approximation. 
One would have to start looking further to the derivatives of A (i.e. F^ u ). In the second 
case, however, scalar background fields give rise to the Higgs mechanism. Except for ad- 
joint color directions that commute with A x , the gauge fields A y and A z will develop large 
masses of order gA x . If A x is large enough, this could plausibly suppress all fluctuations 
of the fields except in commuting color directions. If the unsuppressed fields all commute, 
then non-abelian interactions cannot prevent instability growth. This is equivalent to the 
"abelianization" conjecture of Ref. 0, which supposed that, once non-abelian interactions 
became important, their dynamics would approximately abelianize the fields and so allow 
for continued instability growth. What this conjecture did not appreciate was the difference 
between scalar fields, such as arise in the dimensionally reduced theory, and true gauge 
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fields. 

To make this difference more concrete, we will now analyze a particular example of 
possible relevance to the physics of Weibel instabilities. 



B. An Example: Nielsen-Olesen Instabilities 

Three years ago, before any results had been obtained from three-dimensional simulations, 
Berndt Miiller 18j suggested that three-dimensional instability growth would eventually 



have to stop because, even if the gauge fields did abelianize, Nielsen-Olesen instabilities [19 
would eventually destroy nearly-abelian configurations as the fields continued to grow. In 
this section, we will briefly review these instabilities and discuss how they are suppressed if 
some components of gauge fields are replaced by scalar fields due to dimensional reduction. 

For definiteness and simplicity, let us focus on SU(2) gauge theory. Suppose that there 
is a large magnetic field whose color lies in an Abelian U(l) subgroup. In our application, 
imagine that this is a large, abelianized magnetic field that might have formed due to the 
instability and non-abelian interactions. Over small enough distance and time scales (to be 
discussed later), we can approximate this magnetic field as constant. So, for instance, 

B?(x,t) = B 5 a3 . (2.2) 

Now consider fluctuations about this background field. If we ignore the effects of the hard 
particles represented by Sf(x,p,t) in our application, but just consider pure Yang-Mills 
theory of the soft fields, then this is the problem studied by Nielsen and Olesen [19j. In 
general, the equation of motion for fluctuations a about a classical background A is 

[D V - 2igF^] ab a b u = (2.3) 

if one works in background gauge (D^a^ = 0) and if A satisfies the Yang-Mills equations 
of motion [which is the case for (12.21) ]. Here D = d — igA is the background covariant 
derivative in adjoint color representation. For the background f)2.2p . this equation can be 
decomposed into different sectors classified by (i) the color charge Q = or ±1 of a under 
the color generator T 3 and (ii) its spin projection m s = or ±1 along the Bq axis. Eq. (12. 3p 
becomes 4 

{D 2 + 2m s QgB ) a Qms = (2.4) 

where D = d — iQgA is an Abelian background derivative with charge Qg. For Aq = 0, 
Fourier transforming t to u gives 

[-(V - iQgA) 2 - 2m s QgB ] a Qm = u 2 a Qms . (2.5) 

Mathematically, this is a relativistic generalization of the Schrodinger equation for a particle 
with spin in a constant magnetic field. It has precisely the form of the non-relativistic 
Schrodinger equation with the non-relativistic energy E replaced by u 2 /2M . The 2m s QgB 
term represents the effect on the energy of the interaction between the gluon spin and 



Specifically, a al (x) = J^Q.m, a Q,m a (%) Cq&,„) where Q and m s run over and ±1; C±i = £±i = (1, ±i,0), 
Co =£,o = (0,0,1); and a_ Q _ m . s = a* Q 
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the magnetic field. The values of uo 2 are determined by (12.51) to be the eigenvalues of the 
Hamiltonian 

H Q , ms = (p - QgA) 2 - 2m s QgB . (2.6) 

For the three-dimensional case, one may now follow the discussion of Landau levels in any 
quantum mechanics textbook. The eigenvalues are 

= p\ + (n + \)2gB - 2m s QgB . (2.7) 

The first term is from the component p\\ of momentum parallel to the magnetic field Bq, 
the second is the Landau level energy in the plane orthogonal to B Q , and the last is the 
spin interaction with B . The interesting case here is Q — m s — ±1, for which the lowest 
Landau orbital gives 

u 2 = p\ - gB , (2.8) 

which is negative for pu < (gB Y^ 2 . There are therefore unstable modes in the relativistic, 
classical field theory problem. 5 This is the Nielsen-Olesen instability. Since it produces 
growth of gauge field fluctuations a with non-trivial T 3 charge — that is, of colors which 
do not commute with the original field A — it is an instability which will tend to destroy 
abelianization of the gauge fields. 

But now let's consider the same Hamiltonian (12. 6p for the dimensionally- reduced theory. 
For example, p x = in a two dimensional theory of the yz plane, and so (focusing on the 
previously unstable case Q = m s = ±1): 

H±i,±i = (gA x ) 2 + (p y t 9 A y ) 2 + ( Pz T gA z f - 2gB . (2.9) 

This Hamiltonian is bounded below by (gA x ) 2 — 2gB . If the A x field is large enough that 

(gA x ) 2 > 2gB , (2.10) 

then there will be no negative eigenvalue and so no Nielsen-Olesen instability. In fact, the 
(gA x ) 2 contribution to the relation that determines u 2 is just the Higgs mechanism due to 
a background scalar field A x , as discussed earlier. 
The one-dimensional case is similar, with now 

#±i,±i = H) 2 + (gA y ) 2 + ( Pz T gA z f - 2gB (2.11) 

bounded below by (gA x ) 2 + (gA y ) 2 - 2gB . 

We can usefully recast some of our conditions for the application to Weibel instability 
growth. By treating the magnetic field as constant, we found Nielsen-Olesen instabilities 
in three dimensions with typical moment p» ~ (gB) 1 ^ 2 and growth rates T = 
(gB) 1 / 2 . The constant magnetic field approximation will be okay if these scales are large 
compared to the typical momentum scale k of B. So we need (gB) 1 ^ 2 3> k, which is 

k 2 

fi>— . (2.12) 

g 



5 Note that there is nothing unstable in the non-relativistic Schrodinger problem. The instability uo' 2 < 
in the relativistic, classical held theory problem simply corresponds to a negative energy state in the 
corresponding non-relativistic Schrodinger problem. 
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This is simply the condition for the field to have grown past the point where non-Abelian 
interactions become important. 

For the two-dimensional theory, we can recast condition (12.101) for suppressing Nielsen- 
Olesen instabilities into a similar form. Using B y = D Z A X ~ kA x , the condition is 



So, unless B y is an unusually small component of the magnetic field (because we set A x = 0, 
for example), then the Nielsen-Olesen instability will be suppressed in the dimensionally 
reduced theory just when it first might have become important, as determined by (I2.12p . 

III. CONCLUSIONS 

It is sometimes desirable to simplify a theory, such as through dimensional reduction, 
in order to reduce the cost of numerical simulations. For studying the fate of non-abelian 
Weibel instabilities, we have found qualitatively different behavior in lower- dimensional sim- 
ulations when polarizations A of the gauge field are included that are not gauge fields in 
the lower-dimensional theory. This rules out using one- dimensional (1D+3V) simulations to 
investigate the fate of non-abelian Weibel instabilities. However, two-dimensional (2D+3V) 
simulations can be made to behave similar to the three-dimensional case if the out-of-plane 
polarization of A is discarded. So far, we have not investigated the case of extremely 
anisotropic distributions. 
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APPENDIX A: SIMULATION DETAILS 
1. Basics 

We use the same simulation method as Ref. [1], discretizing dependence of W(x,v,t) 
on velocity v in terms of spherical harmonics Yi m (v) for m < I < / max - For our initial 
conditions, we choose a seed magnetic field configuration and start with zero electric field 
and hard particle fluctuation W. For simulations with small initial conditions, the initial 
magnetic field was chosen such that the gauge fields A have a Gaussian distribution with an 
exponential fall-off in k as in Ref. (8[ Our canonical choice of parameters for two-dimensional 
simulations was lattice spacing a = 0.125/moo on a 256 x 256 lattice of physical size L 2 = 
(256a) 2 = (32/moo) 2 , and Z max = 12 with damping on large / as in Ref. [ll| . 




A relatively simple way to simulate the out-of-plane polarization A x in our 2D+3V sim- 
ulations is to simply run 3D+3V simulations on 1 x N y x N z lattices. The links U x in the 
short periodic direction of size N x = 1 encode the physics of A x (y, z, t). However, we found 




(2.13) 
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it simpler to adapt three-dimensional simulations by instead simulating 2 X N y x N z lattices 
with initial conditions that are translation invariant in the short direction N x = 2. The 
reason is that the procedure for updating the W field, as described in Sec. V B and V D of 
Ref. [201, involves alternating updating W on only the even spatial sub-lattice on one time 



step and only the odd spatial sub-lattice the next time step. For each update of W at a site, 
the algorithm assumes that W at neighboring sites is fixed. That assumption would fail for 
a 1 X N y X N z lattice because every site is its own neighbor. 6 



Initial conditions 



The simulations discussed in the main text all started with tiny seeds for the unstable 
modes. They all have initial exponential growth, as predicted by a perturbative analysis 
of the instability. In Ref. [11], it was noted that linear energy growth sets in immediately 
in three-dimensional simulations if one instead starts with non-perturbatively large seeds 
for the unstable modes. To check the robustness of our conclusions, we have made similar 
simulations in two dimensions. The results are shown in Fig. [5] and [6j As with the case of 
small initial conditions, the late-time behavior is exponential if A x is included and linear if 
not. 




A = A^ty, z, t), M- = x,y,z,t 
A=AJy, z, t), fi = y,z,t, A s = 
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FIG. 5: As Fig. [2] but for two-dimensional simulations with large initial conditions. The initial 
condition is selected from a thermal ensemble with temperature T = 2m 00 /5 2 : , which is then gauge- 
invariantly smeared to remove high-momentum components, exactly as in Ref. [11]. [As in Ref. 
Ill ] , the amount of smearing was chosen so that perturbatively it would correspond to multiplying 
the thermal spectrum for A(k) by exp(— /c 2 /4m^ ).] For the A x = simulation, the links U in the 
x direction are then set to unity. 



Ignoring this issue and proceeding with a periodic lattice with N x = 1 introduces artificial numerical 
instabilities unless one modifies the W update algorithm. This numerical instability can be understood 
analytically by considering a free W field (d t W + v ■ = 0), finding the corresponding discrete-time 
lattice dispersion relation for a given algorithm, and checking whether the solution for the frequency ui 
ever produces a positive imaginary part. 
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FIG. 6: As Fig. 0] but the initial conditions for the 2-dimensional gauge fields A y and A z are large, 
while the out-of-plane component A x is initialized small. 



3. Systematic errors 

Ref. jl] provides many checks that their three-dimensional simulation results are not 
significantly affected by discretization effects. Because the canonical simulation parameters 
we use for two-dimensional simulations in this paper are as good or better than the canonical 
parameters of Ref. one might expect that there are no significant discretization effects. To 
be thorough, however, we present studies of lattice spacing dependence, volume dependence, 
and Z max dependence in Figs. ITHTOl We show here results for simulations with large, non- 
perturbative initial conditions, such as discussed above; results are similar for small initial 
conditions. 
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FIG. 7: Magnetic energy vs. time showing results for several different choices of lattice spacing for 
a volume L 2 = (32/moo) 2 with all polarizations retained (i.e. including A x ). 

The curves in Fig. [7] are like the solid line in left-hand plot of Fig. [3J where we saw 
continued exponential growth, except that here we show explicitly what happens when the 
fields get so large that they become limited by the effects of the discrete lattice. The proof 
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FIG. 8: Linear plot of Magnetic energy vs. time showing results for several different choices of 
lattice spacing for a volume L 2 = (32 /m^) 2 with out-of-plane polarization A x not included. 



that this limit is a lattice artifact is that it increase steadily as the lattice spacing is reduced, 
as seen in Fig. [7J In contrast, the curves in Fig. [S]show the case A x = 0. These curves are 
like the dashed line in the right-hand plot of Fig. where we saw late-time linear energy 
growth. We see here that this behavior is robust as we decrease the lattice spacing, and that 
a = 0.25/moo is an adequately small spacing for our simulations. 

Fig. [9] shows an example of dependence on physical two-dimensional volume L 2 for the 
case of continued exponential growth. (The apparent end to exponential growth shown in 
this figure is the same finite-spacing lattice artifact seen in Fig. [3 We see that nothing 
changes qualitatively as the volume is increased: the exponential growth beyond the non- 
perturbative scale remains. In examining the quantitative variation between the curves, one 
should keep in mind that it is not possible to use the same initial conditions for simulations 
with different volume: in addition to systematic effects, there is statistical variation reflecting 
the random choice of initial conditions. 
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FIG. 9: Magnetic energy vs. time showing results for several different physical volumes with lattice 
spacing a = 0.25/moo . A x is included in these simulations. 

Finally, Fig. [10] is the analogous figure for Z max dependence. We see that Z max = 12 is an 
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adequate approximation to the large Z max limit for our simulations. 
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FIG. 10: Magnetic energy vs. time showing results for several different choices of l ma x with a = 
0.25/moo and I? = (32/rrioo) 2 . A x is included in these simulations. 
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